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Introduction 


The  movement  of  Arctic  ice  floes  is  receiving  greater  interest  today  than 
ever  before.  While  the  tracking  of  ice  floes  is  important  to  aiding  in  Arctic 
navigation,  a  knowledge  of  the  dynamics  of  sea  ice  is  also  vital  to  understanding 
critical  environmental  issues  such  as  global  warming. 

Although  for  many  years,  ice  observations  were  obtained  primarily  under 
proprietary  government  studies,  several  marginal  ice  zone  experiments  have 
recently  been  conducted  which  have  provided  accessible,  enlightening  data  of  ice 
floe  tracks,  ice  compositions  and  floe  geometries.  In  addition  to  large  projects, 
such  as  the  Marginal  Ice  Zone  Experiment  (MIZEX),  the  Coordinated  Eastern 
Arctic  Experiment  (CEAREX),  and  the  Antarctic  Marine  Ecosystem  Research  in 
the  Ice  Edge  Zone  (AMERIEZ),  there  have  also  been  many  smaller  studies.  The 
relatively  new  data  from  these  studies  has  made  it  possible  to  more  accurately 
model  ice  floe  movement  physically  and  numerically.  The  effect  is  actually 
two-fold,  since  the  site  data  has  provided  critical  information  needed  to  prepare 
and  perform  relevant  model  tank  tests.  The  availability  of  the  tank  test  data,  in 
addition  to  the  site  data,  allows  validation  and  further  refinement  of  numerical 
prediction  methods. 

One  of  the  primary  difficulties  in  modelling  (both  physically  and 
numerically)  ice  floe  movement  is  the  numerous  independent  parameters  which 
significantly  affect  floe  behavior.  In  addition  to  geometric  parameters,  such  as 
length,  beam,  draft,  and  under-ice  topography,  the  salinity  profile  of  the  water 
body  in  which  the  ice  floe  moves  has  been  shown  to  have  a  large  effect  on  the 
ice-water  resistance,  and  thus  the  floe’s  motion  (e.g.  Morison,  et.  al.,  1987). 

The  tests  described  in  this  thesis  investigate  the  additional  resistance  that 
internal  wave  creation  induces  on  an  ice  floe  moving  in  a  two-layer,  density 
stratified  water  body.  Stratified  conditions  exist  in  what  is  referred  to  as  the 
Marginal  Ice  Zone  (MIZ)  during  the  spring  melt  season.  The  stratification  is 
created  when  the  ice  floes,  which  have  very  little  salt  content,  melt  and  create  a 
layer  of  relatively  warm,  low  salinity  (and  low  density)  water  on  top  of  the 
existing  water,  which  has  typical  open-ocean  salinity,  Arctic  temperature,  and 
thus  higher  density  values.  The  stratification  is  remarkably  stable,  with  a  sharp 
pycnocline  lasting  for  long  periods  of  time. 


5 


Previous  Laboratory  and  Field  work 


Laboratory  Experiments 

There  have  been  a  limited  number  of  laboratory  experiments  undertaken 
which  investigate  the  added  resistance  due  to  internal  wave  creation  on  moving 
bodies.  However,  each  experiment  described  here  contains  useful  information  for 
further  studies  which  incorporate  similar  tank  experiments. 

Perhaps  the  most  well-known  laboratory  experiments  of  the  type  were 
conducted  by  Ekman  (1904).  His  study  was  motivated  by  field  observations  of  a 
phenomenon  that  he  referred  to  as  “dead  water.”  The  study  illustrated  and 
explained  the  vertical  momentum  transfer  from  a  body  moving  on  the  surface  to 
the  generation  of  internal  waves.  It  was  not  until  much  later  that  similar  studies 
were  conducted  and  their  results  published. 

Muench  and  Hachmeister’s  laboratory  study  (1984)  modelled  sea  ice  keels 
in  a  stratified  tank,  15  m  long,  0.5  m  wide,  and  1  m  deep.  The  study  was  similar 
to  the  one  discussed  in  this  report;  the  additional  drag  that  internal  wave 
generation  induced  on  the  travelling  ice  keel  was  extracted  and  plotted  versus 
non-dimensional  velocity,  U/C,,  where  C,  is  the  first-mode  interfacial  phase 
speed.  (This  is  discussed  later  in  the  Analysis  section.)  Some  of  the  significant 
differences  between  Muench  and  Hachmeister’s  study  and  the  Davidson 
Laboratory  study  include  the  modelled  floe  geometries,  as  well  as  the  parameters 
varied  during  the  tests. 

Ma  and  Tulin’s  (1992)  study  investigated  the  internal  wave  wake 
geometry  and  amplitude  generated  by  a  spheroid  (simulating  a  ship)  moving 
faster  than  the  fastest  internal  wave  phase  speed  (the  supersonic  case)  in  a 
stratified  tank,  12  ft  (3.7  m)  long,  8  ft  (2.4  m)  wide,  with  a  2  ft  (0.61  m)  depth. 
Although  the  relevance  of  this  study  to  Arctic  ice  floe  movement  is  limited,  the 
report  contains  some  interesting  experimental  setup  information. 

Helfrich  and  Melville’s  (1985)  paper  discusses  an  experimental  and 
theoretical  study  of  long,  nonlinear  internal  waves  over  slope-shelf  topography. 
The  laboratory  experiments  were  conducted  in  a  salt-stratified  tank,  24  m  long, 
0.38  m  wide  with  a  0.6  m  water  depth.  Although  the  waves  in  this  experiment 
were  created  by  an  internal-wave  maker,  not  a  moving  disturbance  on  the  free 
surface,  the  information  contained  in  the  report  may  be  quite  valuable  for  further 
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studies.  The  shadowgraph  photographs  of  the  internal  waves  are  particularly 
impressive.  The  experimental  results  compared  very  well  with  their  KdV  model 
predictions.  However,  problems  with  use  of  the  theoretical  model  to  study  actual 
oceanographic  situations  are  discussed. 


Field  Studies 

As  stated  earlier,  there  now  exist  a  great  deal  of  field  observations  which 
may  be  applied  toward  modelling  ice  floe  dynamics.  The  studies  mentioned  here 
represent  a  mere  sampling  of  the  literature  that  has  been  published  to  date. 

The  study  by  Morison  et.  al.  ( 1987)  is  perhaps  one  of  the  most 
comprehensive  studies  of  the  parameters  affecting  ice  floe  drag  in  the  summer 
MIZ  of  the  Greenland  Sea.  The  study,  performed  during  MIZEX  ‘84,  found  the 
ocean  conditions  to  be  highly  variable.  Data  collected  from  the  two  drifting  ice 
stations  include  observations  of  ocean  boundary  layer  turbulence,  profiles  of 
water  temperature,  salinity,  and  velocity,  as  well  as  ice  ablation,  ice 
concentration,  solar  radiation  and  spectral  albedos.  The  stratification’s  impact  on 
the  ice  floe  drag  and  heat  transfer  is  discussed.  It  was  found  that  for  high 
stratification,  the  drag  increased  dramatically,  possibly  due  to  internal  wave 
generation  by  ice  bottom  roughness.  Correlation  with  Muench  and  Hachmeister’s 
(1984)  results  is  also  discussed.  In  the  report  it  is  stressed  that  more  laboratory 
experiments  need  to  be  conducted  to  provide  more  data  to  validate  theoretical 
models.  McPhee  and  Kantha’s  (1989)  paper  later  discussed  these  field 
observations  in  their  study  of  internal  wave  generation  by  sea  ice. 

Sandven  and  Johannessen’s  (1987)  paper  on  internal  wave  observations  in 
the  MIZ  states  two  possible  internal  wave  generators:  tidal  oscillation  over  a 
topographic  feature,  or  ice  keels  moving  in  relation  to  the  water.  The  internal 
wave  characteristics  and  wave  energy  levels  of  the  observed  waves  are  discussed; 
and  a  dispersion  relation  is  solved  numerically  and  compared  with  the 
observations.  They  conclude  that  the  internal  waves  observed  may  be  modelled, 
to  a  good  first  approximation,  as  horizontally  propagating  plane  gravity  waves  in 
an  inviscid  fluid. 

There  are  also  several  field  studies  of  internal  wave  dynamics  outside  of 
the  Marginal  Ice  Zone.  For  example,  Watson,  et.  al.  (1992)  discuss  measurements 
of  a  novel  experiment  in  which  three  ships  were  used  to  generate  internal  wave 
wakes  in  a  stratified  sea  loch.  Valuable  results  are  obtained,  particularly  in 
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reference  to  the  dependence  of  disturbing-body  length  to  generated  internal  wave 
characteristics. 
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Previous  Theoretical  and  Numerical  Work 


Internal  Wave  Theory 

The  effect  of  internal  waves  on  ice  floe  movement  is  discussed  by  Bruno 
(1991).  The  assumptions  in  his  development  of  the  internal  wave  equations  are 
applicable  to  the  study  addressed  in  this  thesis.  Therefore,  the  rigorous  derivations 
provided  of  the  equations  governing  internal  wave  generation  should  provide 
guidance  in  further  analyses  of  the  results  found  in  the  Davidson  Laboratory 
study. 


Yih  (1960)  presents  a  model  of  wave  motion  in  a  stratified,  inviscid  fluid 
obtained  by  reducing  the  governing  equations  to  a  Sturm-Liouville  system.  The 
paper  provides  applications  for  continuously  stratified  fluids,  as  well  as  those 
containing  density  discontinuities.  The  result  is  a  very  versatile  numerical  model 
with  broad  applicability. 

Zhu,  et.  al.  (1987)  present  and  discuss  two  theoretical  models  of 
dispersive  long  waves  propagating  in  a  stratified  fluid.  One  is  an  inhomogeneous 
KdV  (IKdV)  model,  while  the  other  is  a  Boussinesq-class  model.  Both  models 
were  developed  for  predicting  the  internal  waves  generated  by  a  pressure 
disturbance  moving  in  incompressible,  inviscid  and  stratified  fluid  systems.  The 
two  methods  are  compared  to  each  other  and  to  experimental  data.  Although  their 
discussions  of  two-layer  stratified  systems  is  limited,  their  findings  show  potential 
for  application  to  this  ice  floe  study. 

There  are  also  many  other  studies  which  investigate  internal  wave  energy 
and  internal  wave  creation  by  specific  types  of  forcings.  Ball’s  (1963)  paper 
discusses  the  transfer  of  a  surface  wave’s  energy  to  an  internal  wave,  with 
detailed  derivations  for  the  shallow  water  case.  It  may  be  possible  to  draw  an 
analogy  of  the  surface  wave’s  impartation  of  energy  to  an  internal  wave  field  to 
that  of  an  ice  floe.  However,  complications  arising  due  to  wave  field  coupling, 
which  are  nonexistent  in  ice  floe  dynamics,  may  hinder  comparisons  of  the  two 
phenomena. 
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Ice  Floe  Numerical  Modelling 


Several  ice-ocean  numerical  models  have  been  developed.  The  types  of 
models  vary  widely,  from  short-term  (order  of  hours  or  days)  to  long  term  (order 
of  months  or  years).  However,  none  have  included  the  additional  resistive  and 
coupling  effects  attributable  to  internal  wave  creation. 

Some  of  the  recent  models  which  incorporate  ice-ocean  coupling  effects 
include  Bruno  and  Madsen  (1989),  Kantha  and  Mellor  (1989),  and  Pritchard  et. 
al.  (1990).  Coupling  of  the  ice  and  ocean  affects  the  ice  motion  and  the  dynamics 
of  the  upper  ocean.  The  ice-water  resistance,  which  is  the  source  of  the  coupling, 
remains  largely  unresolved.  The  resistance  at  the  ice-water  interface  is  generally 
attributed  to  only  two  factors:  skin  friction  resistance  associated  with  the 
small-scale  roughness  of  the  under-ice  topography,  and  form  drag  associated  with 
pressure  disturbances  from  large-scale  features  such  as  keels.  However,  in  a 
stratified  ocean,  another  coupling  mechanism  due  to  the  generation  of  internal 
waves  exists.  It  has  been  shown  that  this  phenomenon  may  significantly  affect  the 
ice-water  resistance  and  thus  the  ice  motion  (e.g.  Muench  and  Hachmeister, 

1984). 


Separate  studies  of  ice-water  drag  coefficient  determination  for 
implementation  into  ice  floe  numerical  models  also  neglect  internal  wave 
generation  effects.  From  studies  such  as  Madsen  and  Bruno’s  (1987)  relatively 
simple  methodology  for  calculating  ice  floe  drag  coefficients  to  McPhee’s  (1982) 
thorough  analysis  of  ice-water  dynamics,  explicit  incorporation  of  internal  wave 
generation  is  not  identified.  In  McPhee’s  (1982)  report,  a  rigorous  analysis  of 
several  parameters  affecting  ice-water  drag,  including  boundary  layer  alteration 
due  to  rapid  melting,  is  addressed.  McPhee  discusses  deficiencies  in  ice-water 
drag  treatment  of  other  numerical  models,  and  provides  several  methods  for 
incorporating  the  momentum  flux  between  the  ice  and  the  water  in  those  models. 
However,  later  studies  indicate  the  necessity  of  also  including  internal  wave 
generation  effects  (e.g.  McPhee  and  Kantha,  1989). 

Although  discussion  of  the  numerical  models  in  depth  would  transcend  the 
scope  of  this  thesis,  it  is  important  to  emphasize  that  internal  wave  effects  on 
ice-water  drag,  and  therefore  predictive  models,  still  demand  further  attention. 
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The  Davidson  Laboratory  Project 


Introduction 

The  series  of  experiments  discussed  in  this  thesis  investigated  the 
increased  drag  exhibited  by  ice  floes  in  salinity-stratified  waters.  The  tests 
examined  the  dependence  of  the  additional  ice-water  drag  and  the  characteristics 
of  the  generated  internal  wave  on  parameters  such  as  floe  velocity,  floe  draft, 
upper  layer  thickness,  and  under-ice  topography.  These  parameters  were  therefore 
varied  in  the  test  matrix  to  achieve  information  on  the  corresponding  trends.  Both 
non-strati  fled  and  stratified  tests  were  performed. 

The  tests  were  conducted  in  Davidson  Laboratory’s  Tank  3  Facility.  The 
tank  is  313  ft  (95.4  m)  long,  12  ft  (3.6  m)  wide,  and  permits  water  depths  up  to  6 
ft  (1.8  m).  The  tank  was  designed  primarily  as  a  high-speed  naval  architectural 
towing  tank.  However,  many  other  hydrodynamic  projects,  including  offshore 
structure  wave  impact  studies,  reef  flow  visualizations  and  coastal  beach  erosion 
studies  have  been  conducted  in  this  tank.  The  towing  carriage  permits  velocities 
up  to  100  ft/s  (30  m/s)  with  a  0.01  ft/s  (0.003  m/s)  resolution.  The  towing 
mechanism  provides  extremely  accurate  and  consistent  low  speed  towing 
velocities. 

For  the  stratified  experiments,  the  tank  was  salinity-stratified  by  creating  a 
saline  water  body  and  layering  fresh  water  above  the  salt  water.  The  stratification 
desired  was  a  two-layer  configuration,  to  simulate  marginal  ice  zone  conditions 
observed  by  Morison  et.  al.  The  stratification  procedure  is  described  and  a  sample 
achieved  salinity  profile  is  provided  in  the  Large-Scale  Experimental  Setup 
section. 


The  ice  floe  models  were  constructed  of  extruded  polystyrene  and  fitted 
with  an  interior  wood  frame  to  facilitate  connection  to  the  towing  apparatus.  A 
stiff  spring  balance  with  a  maximum  capacity  of  10  lb,  used  to  obtain  drag 
measurements,  connected  the  floe  to  the  towing  apparatus.  Internal  wave  data  was 
acquired  through  two  methods:  time-history  wave  wire  records  and  VHS  video 
tape  recordings.  (Both  are  also  described  further  in  the  Large-Scale  Experimental 
Setup  and  Instrumentation  subsection.)  See  Figure  1  for  sketch  of  the  test  setup. 


II 


Scaling 

Standard  naval  architectural  practices  were  followed  in  the  scaling  of  the 
ice  floe  tests.  Froude  number  similarity  between  full  scale  and  model  scale  floes 
was  assumed.  Bruno  (1992)  discusses  the  impracticalities  of  Reynolds  scaling, 
and  gives  further  justification  for  employing  Froude’s  hypothesis  to  this  series  of 
experiments. 

The  full  scale  floes  were  assumed  to  have  the  following  idealized 
geometric  characteristics: 

length,  L=  100  m 
beam,  B=  1 00  m 
draft,  T=  8.3  to  16.6  m. 

A  length  scale  ratio  of  109  was  used,  which  yielded  model  floe  dimensions  of: 

length,  L=  3  ft  (0.91  m) 
beam,  B=  3  ft  (0.9 1  m) 
draft,  T=  3  to  6  in  (7.6  to  15.2  cm). 
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The  majority  of  the  upper  layer  thicknesses  modelled  corresponded  to  full  scale 
thicknesses  between  19  m  and  25  m,  with  one  set  of  tests  performed 
corresponding  to  a  full-scale  upper  thickness  of  33  m.  Although  the  full-scale  ice 
drafts  are  high  compared  to  marginal  ice  zone  observations  (Bourke  and 
McLaren,  1992),  preliminary  tank  tests  indicated  that  these  corresponding  higher 
model  drafts  were  required  to  create  quantifiable  drag  forces  and  internal  wave 
disturbances  in  the  experiments.  In  any  case,  it  appears  that  the  critical  parameter 
is  not  the  absolute  floe  draft,  but  rather  the  draft  relative  to  the  upper  layer 
thickness,  which  was  representative  of  full-scale  conditions. 

Non-dimensional  total  resistance  coefficients  were  obtained  by  applying 
the  following  formula: 

C  = 

'  pAlP' 

where  F,  is  the  total  measured  drag  force,  A  is  the  wetted  surface  area,  and  U  is 
the  towing  speed.  (Note  the  variation  from  some  conventions  of  defining 


Pursuant  to  common  naval  architectural  practice,  the  total  drag  force 
acting  on  the  ice  floe  was  separated  into  two  components:  a  frictional  resistance 
component,  Fp  and  a  residual  resistance  component,  F,: 

F,=  Ff^  Fr. 

The  residual  resistance  component  includes  the  effects  from  all  of  the 
non-frictional  drag  contributors.  The  coefficients  and  are 
non-dimensionalized  in  a  similar  manner  as  C,,  employing  the  respective  resistive 
component  forces,  Ff  and  F^  Therefore,  the  total  resistance  coefficient  may  be 
written: 

Ct  —  Cf  +  Cr 


Arctic  ice  floe  observations  have  indicated  that  a  moving  ice  floe  creates 
extremely  little  to  no  measurable  surface  wave  creation.  The  Froude  number  of 
the  model  tests  was  therefore  limited  to  avoid  creation  of  surface  waves  during 
the  tank  experiments.  This  condition  then  simplified  the  ice-water  drag  analyses 
by  allowing  assumption  of  zero  surface  wavemaking  drag. 

In  the  unstratified  tests,  the  only  sources  of  resistance  were  therefore 
assumed  to  be  skin  friction  and  form  drag,  hence  the  value  of  F^  is  equal  to  the 
form  drag,  F^,  i.e.  F_.  =  Fj.  Therefore,  for  unstratified  tests,  the  above  equation 
may  be  written: 


P I.  unstratified  —  Pf  P d  • 
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It  is  recognized  that  most  conventions  non-dimensionalize  the  form  drag 
coefficient,  Cj,  using  the  frontal  projection  area  (i.e.  beam  x  draft).  However,  in 
this  thesis,  continuity  throughout  all  coefficients  is  preserved  by 
non-dimensionalizing  using  the  total  wetted  surface  area.  To  avoid  any 
confusion,  the  coefficient  Cj  will  not  be  used,  and  the  form  drag  component  of 
the  total  resistance  will  hereafter  be  referred  to  only  as  the  “unstratified  residual 
resistance,”  and  the  corresponding  coefficient  written  as  “C  j.”  Therefore: 

C t,  unstratified  ~  Cy  +  C r,  unstratified  • 

In  the  stratified  tests,  the  residual  resistance  included  the  additional  drag 
caused  by  internal  wave  creation,  i.e.  +  F^,,,  where  F;,,  is  the 

internal  wave  drag  component.  Therefore  the  total  drag  coefficient  for  the 
stratified  tests  was  written: 

Cf  stratified  Cfi  +  Cr.  unstratified  Ctw. 


Flat  Plate  Tests 

Before  the  ice  floes  were  tested,  flat  plate  tests  were  performed  to 
determine  the  skin  frictional  resistance  of  the  material  used  to  construct  the  model 
ice  floes.  The  “flat  plate”  used  in  such  a  test  is  actually  a  streamlined,  thin  board 
of  the  material  of  interest  with  tapered  ends,  which  is  mounted  vertically  and 
towed  lengthways  so  as  to  minimize  form  drag  and  surface  wave  production.  The 
assumption  in  flat  plate  tests  is  that  the  resistance  due  to  skin  friction  is  so 
proportionally  greater  than  the  other  components  of  the  resistance  (such  as 
wavemaking  and  form  drag),  that  the  total  resistance  measured  is  entirely 
attributable  to  skin  friction  resistance.  Towing  velocities  are  also  limited  so  as  to 
prevent  surface  wave  creation. 

Following  these  assumptions,  a  flat  plate  was  created  out  of  the  material 
used  to  construct  the  ice  floe  models,  and  attached  to  the  carriage  with  a  draft  of  9 
in  (23  cm).  (A  sketch  of  the  flat  plate  used  for  these  tests  is  shown  in  Figure  2.) 

A  soft  spring  drag  balance  with  a  maximum  capacity  of  10  lb  and  a  resolution  of 
0.01  lb  was  used  to  measure  the  plate’s  drag  force.  The  friction  coefficient  for  the 
floe  material  was  determined  from  the  flat  plate  tests  as: 

^  _  F toial,  flat  plate 

pAlF  ’ 

where  F,^,  was  the  total  measured  drag  force  from  the  flat  plate  test. 
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Figure  2:  Flat  Plate  Model  -  Plan  and  Profile  Views 


The  Cf  values  varied  only  slightly  during  the  experiments,  and  an  average 
value  of  Cf=  4.5  x  10'^  was  adopted  for  the  speed  range  of  0.3  to  0.6  m/s,  which 
corresponds  to  a  Reynolds  number  range  of  3x10^  to  6x10^.  This  skin  friction 
information  was  later  used  to  extract  the  residual  resistance  component  from  the 
total  resistance. 


Small-Scale  Experiments 

A  series  of  relatively  small-scale  qualitative  experiments  were  conducted 
before  the  main  tests  in  a  Plexiglass  tank  (unofficially  referred  to  as  the 
“Davidson  Laboratory  Tank  4  Facility”)  measuring  8  ft  (2.4  m)  long  by  2  ft  (0.61 
m)  wide,  with  a  variable  water  depth  up  to  1.8  ft  (0.55  m).  (See  Figure  3.)  This 
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Figure  3:  Small-Scale  Experimental  Tank 


setup  was  ideal  for  experimenting  with  different  methods  of  creating  a  stratified 
water  column,  as  well  as  with  methods  for  sharpening  the  interface  between  the 
two  salinity  layers.  The  volume  of  the  tank  was  small  enough  to  easily  drain  and 
refill  several  times  per  day  as  needed,  while  still  representing  a  large  enough  scale 
to  provide  valid  information  on  the  behavior  of  the  stratification  and  experimental 
methods  in  the  larger  tank.  This  smaller  tank  also  provided  a  convenient  setup  in 
which  to  experiment  with  different  internal  wave  observation  techniques. 

During  these  tests,  several  methods  of  creating  a  stably  stratified  water 
column  in  the  laboratory  setting  were  tested.  The  techniques  included  several 
methods  which  introduced  fresh  water  on  top  of  a  tank  of  salt  water,  and  some 
which  introduced  salt  water  at  the  bottom  of  a  tank  of  fresh  water.  The  latter 
methods  were  unsuccessful,  due  to  excessive  mixing  of  the  two  fluids  at  even 
very  slow  salt  water  input  flow  rates.  A  very  gradual  halocline  resulted,  which 
was  unacceptable  for  these  experiments. 

The  process  decided  upon  for  application  in  the  larger  tank  is  actually 
quite  simple,  and  was  easily  implemented  in  the  larger  tank.  The  method  requires 
a  salt  water  layer  to  be  established,  and  then  introduces  fresh  water  on  top 
through  hoses  mounted  perpendicular  to  plywood  boards  which  are  floating  on 
the  surface.  The  fresh  water  sheets  off  of  the  boards  in  a  mostly  horizontal 
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direction,  and  layers  on  top  of  the  fluid  beneath  with  very  little  disturbance  to  or 
mixing  with  the  salt  water  layer.  See  Figure  4  for  schematic. 


Fresh  Water  Input 


Hose  Mounting  Bracket 


Fresh  Water 
Layer 


Salt  Water  Layer 


Figure  4;  Sketch  of  Stratification  Process 


It  was  also  during  these  tests  in  the  smaller  tank  that  the  SIS  (Stratification 
Intensification  System)  instrument  was  developed  and  a  small  prototype  was 
built.  The  instrument  intensifies  the  stratification  of  the  water  column  (  sharpens 
the  halocline)  by  removing  a  horizontal  “strip”  of  mixed  salinity  water.  A 
schematic  of  its  operation  is  shown  in  Figure  5.  Figure  6  shows  a  photograph  of 
the  prototype,  better  known  as  “Little  SIS.” 
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Figure  5:  Schematic  of  SIS  operation 


Figure  6:  “Little  SIS” 

Large-Scale  Experimental  Setup  and  Instrumentation 

At  the  commencement  of  the  tests,  the  tank  water  level  was  lowered  to  the 
desired  depth  of  the  salt  water  layer,  4  ft  ( 1 .2  m).  This  water  was  then  salinated  to 
approximately  35  ppt,  by  adding  999  NaCl,  99.9%  pure  table  salt.  The  water  was 
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mixed  thoroughly  by  mounting  a  “mixing  mesh”  on  the  tank  carriage  and  towing 
it  several  times  throughout  the  length  of  the  tank. 

The  method  developed  in  the  small-scale  tests  for  layering  fresh  water  on 
top  of  the  salt  water  proved  to  be  very  successful  in  the  large  tank.  Five  plywood 
sheets,  each  measuring  6  ft  by  12  ft  (1.8  m  by  3.6  m)  were  distributed  along  the 
length  of  the  tank.  A  freshwater  discharge  hose  was  mounted  perpendicular  to 
each  board,  and  fresh  water  was  introduced  until  a  layer  thickness  between  7  and 
9  in  (18  to  23  cm)  was  achieved.  The  resulting  salinity  profile  exhibited  a  sharp 
halocline,  as  desired.  A  sample  salinity  profile  is  shown  in  Figure  7. 


Salinity  Profile 
before  Run  #8 


Figure  7:  Sample  Achieved  Salinity  Profile 


The  “Big  SIS”  was  constructed  of  expanded  PCP,  and  proved  to  be  a 
successful  method  of  sharpening  the  halocline  in  the  larger  tank.  The  instrument 
measured  1 1.5  ft  wide  (3.5  m),  3  ft  (0.91m)  long,  and  3  in  (7.6  cm)  high  at  the 
inflow  (leading)  edge.  (See  Figure  8.)  The  hoses  which  exited  from  the  trailing 
edge  of  the  wedge  were  connected  to  a  manifold  which  was  then  connected  to  a 
pump  suspended  from  the  towing  carriage.  The  discharge  from  the  pump  was 
then  sent  to  a  large  skimmer  drain.  The  wedge  itself  was  attached  to  the  carriage 
and  was  towed  the  length  of  the  tank,  as  needed  to  remove  a  strip  of  mixed 
salinity  water.  The  tow  speed  was  adjusted  to  balance  the  flow  into  the  wedge 
with  the  flow  out  of  the  system  so  as  to  avoid  disturbances.  Symbolically: 

Qin  ~  Qout  ~  ^ frontal,  projected  ^  U, 
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where  was  the  frontal  projected  area  equal  to  the  maximum  (leading 

edge)  thiclmess  of  the  wedge,  and  U  was  the  carriage  towing  speed.  If  necessary, 
additional  fresh  water  was  then  introduced  on  the  surface  (by  the  method 
previously  described)  to  replenish  the  discharged  volume  of  water. 


Figures:  “Big  SIS” 


Several  methods  of  obtaining  internal  wave  recordings  were  tested.  Some 
of  the  unsuccessful  methods  included  dye-injection  techniques,  shadow 
photography  visualization  (from  Mowbray,  1967),  echo-sounder  recordings,  and 
traditional  wave  wire  methods. 

The  problem  with  the  dye-injection  techniques  arose  from  the  high 
diffusivity  of  the  dyes  in  both  fresh  and  salt  water.  No  dye  was  found  that  would 
either  remain  at  the  interface  of  the  two  layers,  or  stay  primarily  in  one  of  the  two 
layers  for  any  reasonable  time  scale.  Reports  of  other  projects  in  which 
dye-injection  techniques  for  internal  wave  visualization  were  successful 
incorporated  density  stratifications  much  greater  (by  at  least  a  factor  of  10)  than 
in  the  Davidson  Laboratory  experiments. 

The  primary  problem  with  implementation  of  Mowbray’s  (1967) 

Schlieren  and  Shadowgraph  techniques  in  the  Davidson  Laboratory  test  series 
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was  the  requirement  of  “observation  windows”  within  the  test  section,  or 
preferably,  an  entirely  glass-sided  tank.  Since  Tank  3  has  solid  concrete  walls  and 
no  observation  windows,  the  numerous  pieces  of  apparatus  required  to  be 
installed  inside  of  the  tank  to  implement  this  method  would  have  greatly 
disrupted  the  internal  wave  motion,  so  this  method  was  therefore  deemed 
impractical.  Successful  implementations  of  this  technique  may  be  seen  in  the 
results  from  Helfrich  and  Melville’s  (1985)  laboratory  experiments. 

An  echosounder  was  tested  as  well,  and  although  this  technique  is 
promising  for  acquiring  internal  wave  data  in  future  studies,  the  particular 
instrument  tested  for  these  experiments  unfortunately  did  not  have  enough 
sensitivity  to  accurately  register  the  internal  wave  location.  An  echosounder  with 
greater  sensitivity  and  resolution  appears  to  be  an  ideal  instrument  for  this  type  of 
test  because  of  its  non-obtrusiveness. 

Traditional  wave  wires  were  also  impractical  for  use  in  this  study,  because 
of  the  severe  corrosive  nature  of  salt  water.  Virtually  every  wave  wire  instrument 
submerged  into  the  tank  corroded  away  within  minutes. 

However,  two  methods  for  internal  wave  observation  proved  to  be 
successful,  and  provided  accurate  recordings  of  the  internal  wave  created  by  the 
moving  ice  floe.  The  first  method  used  an  internal  wave  wire  probe  designed  and 
built  for  this  project,  and  the  second  method  employed  a  submerged  video  camera 
to  record  observations  of  the  internal  wave. 

The  internal  wave  probe  design  was  essentially  a  modified  conventional 
wave  wire  design.  The  instrument  consisted  of  a  pair  of  platinum  wires  22  in 
(55.9  cm)  long,  spaced  0.5  in  (1.3  cm)  apart,  which  were  encased  in  a  1  in  (2.54 
cm)  diameter  hollow  plastic  tube.  The  tube  had  open  slots  which  allowed  free 
vertical  movement  of  the  fresh  water  /  salt  water  interface.  The  internal  wave 
wire’s  operation  is  analogous  to  a  surface  wave  wire’s.  The  internal  wave 
instrument  was  submerged  below  the  free  surface  with  a  portion  of  the  probe 
penetrating  below  the  salt  water  interface.  The  probe  measured  average 
conductivity  of  the  fluid  body,  and  converted  this  into  a  voltage.  Changes  in  the 
level  of  the  salt  water  produced  an  associated  linear  change  in  conductivity,  and 
thus  voltage  reading.  The  output  from  the  instrument  was  then  fed  into  a  strip 
chart  recorder,  and  also  branched  into  the  Davidson  Laboratory’s  computer  data 
acquisition  equipment.  A  schematic  of  the  internal  wave  wire  electronics  is  shown 
in  Figure  9.  A  photograph  of  this  internal  wave  wire  is  shown  in  Figure  10. 


AC/OC  CONVERTER 


ACTIVE  flLlER 


i50  amplifier 


Figure  10:  Internal  Wave  Wire 


Additional  internal  wave  observations  were  obtained  through  a  submerged 
video  camera,  mounted  on  the  side  of  the  tank.  During  the  small-scale  tests,  the 
lower,  saltier  layer  was  clouded  due  to  dust  particles  in  the  salt  which  was  used  to 
salinate  the  water,  while  the  upper,  fresh  water  layer  was  clear.  However,  in  the 
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larger  scale  tests,  the  salt  used  to  create  the  salt  water  was  extremely  pure,  and  the 
layer  was  clear.  However,  after  the  fresh  water  was  layered  above,  it  was 
observed  that  the  tiny  air  bubbles  in  the  “newer”  layer  made  the  water  appear 
cloudy  when  a  bright  light  was  shone  from  above,  while  the  lower,  saltier  layer 
remained  clear.  This  situation  therefore  facilitated  easy  visualization  of  the 
interface  between  the  two  layers  when  viewed  through  a  camera  mounted  below 
the  surface.  This  technique  was  therefore  implemented,  and  provided  excellent 
internal  wave  recordings. 

Difficulties  were  also  encountered  with  salinity  profile  measurement 
instruments.  Although  several  electrical  conductivity  probes  were  designed  and 
built  during  this  project  to  measure  the  electrical  conductivity  and  thus  salinity 
through  the  water  column,  severe  corrosion  and  non-linearity  problems  were 
encountered  which  required  dismissal  of  these  instruments. 

A  successful,  although  laborious,  method  was  ultimately  employed  to 
obtain  the  salinity  profiles.  “Grab  samples”  were  taken  at  different  discrete  levels 
through  the  water  column;  a  hand-held  refractometer  was  then  used  to  measure 
the  salinity  of  each  sample. 
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The  Test  Matrix 


Four  different  ice  floes  were  tested,  each  with  different  under-ice 
topographies  as  described  below: 

Floe  #1:  flat  bottom 

Floe  #2;  low-frequency  sinusoid  bottom  profile 
Floe  #3:  higher-frequency  sinusoid  bottom  profile 
Floe  #4:  irregular  bottom  topography,  typical  of  deformed  first-year  ice 
observations 

The  lower  frequency  sine  wave  of  Floe  #2’s  bottom  topography  had  a 
wavelength  of  1.5  ft  (0.46  m).  The  higher  frequency  sine  wave  of  Floe  #3’s 
topography  had  a  wavelength  of  1 .0  ft  (0.30  m).  Each  sine  wave  had  an  amplitude 
of  1 .0  in  (2.54  cm).  Floe  #4  had  an  irregular  bottom  modelled  after  Lange  and 
Eicken’s  (1991)  observations  of  class  11  ice.  (See  Figures  1 1  through  14  for 
profiles.)  Each  floe  had  a  square-shaped  plan  view,  3  ft  (0.91  m)  side  length,  and 
a  rectangular  frontal  projection  with  a  height  of  1  ft  (0.30  m). 


Figure  1 1 :  Profile  View  of  Floe  #  1 
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Figure  12:  Profile  View  of  Floe  #2 


Figure  13:  Profile  View  of  Floe  #3 
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Figure  14:  Profile  View  of  Floe  #4 


Non-Stratified  Test  Matrix 

The  unstratified  test  matrix  includes  conditions  similar  to  the  stratified 
tests  shown  in  Table  1  for  each  of  the  four  ice  floe  models.  Unfortunately,  many 
of  these  tests  were  not  completed  at  the  time  this  thesis  was  written.  (Floe  #3  and 
Floe  #4  have  not  been  tested  at  all  in  the  non-stratified  condition.)  However,  all 
of  the  remaining  tests  will  be  performed  during  June  of  1993,  and  a  later  paper 
will  present  the  complete  results. 


Stratified  Test  Matrix 

The  stratified  tests  proved  to  be  massively  time  consuming,  due  to  the 
lengthy  stratification  establishment  time  requirements,  and  then  limitations  on  the 
time  during  which  the  tank  remained  suitably  stratified.  Sharpening  the  interface 
through  use  of  the  SIS  was  also  time  consuming.  In  addition,  after  each  run,  very 
long  internal  waves  propagated  along  the  length  of  the  tank;  therefore  more  than  a 
half  hour  time  period  between  the  end  of  one  run  and  the  beginning  of  the  next 
was  required  to  allow  the  internal  waves  to  dissipate. 
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To  simplify  and  streamline  the  test  matrix,  modifications  were  made  to  the 
originally  intended  stratified  test  matrix  during  the  experiment  based  on 
observations  of  the  internal  wave  behavior  and  drag  trends.  The  resulting  test 
matrix  is  displayed  in  Table  1,  with  the  corresponding  run  numbers  for  each  test 
point  identified. 


Table  1:  Stratified  Test  Matrix 
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Test  Results 


Data  Acquisition 

Resistance  and  wave  wire  data  were  acquired  by  the  Davidson  Laboratory 
DAPS  data  acquisition  program.  Data  was  acquired  within  a  50  ft  (15  m)  test 
length  of  the  ice  floe’s  travel.  Liberal  acceleration  and  stabilization  allowances 
were  incorporated  before  this  data  acquisition  length.  The  resistance  values 
obtained  were  calculated  by  averaging  the  input  over  the  test  length.  Time  history 
files  of  the  wave  wire  data  were  also  obtained.  A  summary  of  the  test  data  is 
included  in  Appendix  A.  Strip  chart  recordings  of  the  wave  wire  readings  were 
also  obtained  and  are  included  in  Appendix  B. 

Salinity  profiles  within  the  data  acquisition  portion  of  the  tank  were  taken 
before  and  after  each  test  run.  These  results  are  presented  in  Appendix  C. 

VHS  videotape  recordings  of  the  internal  wave  created  by  the  ice  floe 
were  obtained  by  mounting  a  small  camera  on  the  side  of  the  tank  at  the  level  of 
the  halocline.  Unfortunately,  no  still  photographs  were  taken  during  the  tests, 
since  the  number  of  intrusions  in  the  water  column  was  kept  as  small  as  possible. 


Non-Stratified  Test  Results 

At  this  time,  only  Floe  #1  and  Floe  #2  have  been  tested  in  the 
non-stratified  condition.  Tables  2  and  3  include  the  total  resistance  coefficients 
(^t,uiwraufied  values)  from  Floe  #1  and  Floe  #2,  while  Tables  4  and  5  list  the  residual 
resistance  coefficients  (C,  values)  for  the  two  floes. 
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Table  2:  Unstratified  Test  Total  Resistance  Coefficients,  C,  values, 

from  Floe  #1 


Table  3:  Unstratified  Test  Total  Resistance  Coefficients,  C,  values, 

from  Floe  #2 


Table  4:  Unstratified  Test  Residual  Resistance  Coefficients, 
_  C,  values,  from  Floe  #  1 
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Table  5:  Unstratified  Test  Residual  Resistance  Coefficients 
_ C - -  values,  from  Floe  #2 _ 


Non- Strati  tied  Test 

Residual  Resistance  Coefticienis.  Floe  *f2 


I  Ic*  Floe  I 

Velocity  i 

Draft 

,T 

0  1  tt/s 

0  2fl/s  I 

0  3  ft/S  1 

0  4  ft/S  j 

0  5  ft/s  j 

0  6  ft/s 

0.7  ft/s 

0  8  ft/S  (1 

in 

cm 

3  cm/s 

6  cm/s  i 

9  cm/s  1 

12  cm/s  i 

1 5  cm/s  i 

1 0  cm/s 

21  cm/s 

24  cm/s  1 

2.0 

5  1 

--•-1 

'r'r.-r'4 

■BHII 

— 

3  0 

0.029  1 

0  030  1 

0.029  1 

0.029 

jjssr.- ' 

4  0 

BHI 

... 

mBii 

Bimi 

illlMliiii 

5  0 

0.049  1 

0.044  1 

0.046  1 

0.046 

L60 

15  2 

0.055  i 

0.056  i 

0.056  ! 

0,054 

1  » 

Stratified  Test  Results 

Tables  6  through  9  present  the  total  resistance  coefficients  (C, 
values)  from  the  stratified  tests.  Plots  of  the  data  appear  in  the  Discussion  and 
Preliminary  Analysis  section. 


Table  6:  Stratified  Test  Total  Resistance  Coefficients,  C^  sa^,i^Kd  values, 

from  Floe  #1 


Stratified  Test 

Total  Resistance  Coefficients.  Floe  #1 


Icefloe 

Draft,  T 

in  cm 

T/h, 

Velocity  | 

0  1  tt/s 

3  cm/s 

0.2  ft/s 

6  cm/s 

0,3  ft/s 

9  cm/s 

0  4  h/s 

12  cm/s 

0.5  ft/s 

15  cm/s 

0.6  ft/s 

18  cm/s 

0.7  tt/s 

21  cm/s 

3.0 

7.6 

0.50 

0.054 

0.072 

0.030 

5.0 

12.7 

0.83 

0.090 

5.0 

12.7 

0  67 

u 

0.089 

6.0 

15.2 

0.80 

0.095 

0.140 

6.0 

15.2 

0  67 

0.241 

0.112  1  0.103 
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Table  7:  Stratified  Test  Total  Resistance  Coefficients,  C,  values, 

from  Floe  #2 


Stratified  Test 

Total  Resistance  Coefficients,  Floe  #2 


lc«FkM 

Draft,  T 

in  cm 

T/h. 

Velocity  I 

0.2  ff/s 

6  cm/s 

0.3  tt/s 

9  cm/s 

0.4  tt/s 

12  cm/s 

0.5  tt/s 

1 5  cm/s 

0.6  ft/s 

18  cm/s 

0,7  ft/s 

21  cm/s 

0.8  ft/s  1 
24  cm/s  1 

2.0  S.l 

0  36 

0.027  I  0.020 

0.019  I 

[EEm 

0.55 

0.066  ' 

0.076  1  0,044  i  0.032 

Table  8:  Stratified  Test  Total  Resistance  Coefficients,  C,  values, 

from  Floe  #3 


Stratified  Test 

Total  Resistance  Coefficients.  Floe  #3 


lc«Flo« 

Velocity  j 

Draft 

T/h. 

0.1  ft/s 

0.2  ft/s 

0.3  tt/s 

0  4  ft/s 

0.5  ft/s 

0.6  ft/s 

0.7  ft/s  1 

lin 

cm 

3  cm/s 

6  cm/s 

9  cm/s 

12  cm/s 

15  cm/s 

18  cm/s 

21  cm/s 

EQ 

■0 

0.46 

0.076 

0.052 

0.029 

3  0 

7.6 

0,43 

0.079 

4.0 

10.2 

0.62 

0.202 

0.103 

0.069 

.0.037 

4,0 

10.2 

0,57 

0.113 

Table  9:  Stratified  Test  Total  Resistance  Coefficients,  C,  5^,^  values, 

from  Floe  #4 

Stratified  Test 

Total  Resistance  Coefficients.  Floe  #4 


le«  Fto« 

Velocity  | 

Draft 

T/h. 

0.1  ft/s 

0.2  ft/s 

0.3  ft/s 

0  4  ft/s 

0.5  tt/s 

0.6  ft/s 

0.7  ft/s  j 

I  in 

cm 

3  cm/s 

6  cm/s 

9  cm/s 

■2  cm/s 

■  5  cm/s 

1 0  cm/s 

21  cm/s  1 

^£1 

5.1 

0.29 

1 

0,033 

li 

m 

0,43 

0.044 

■ _ 

1  4.0  10.2 

0.57 

0.060 

0.103  I 

I  0.063  I 

j  0.030  1 

Limited  wave  wire  data  was  recorded.  During  a  few  of  the  tests, 
difficulties  with  the  instrument  were  encountered  due  to  corrosion  of  the  wires. 
However,  valid  results  were  obtained  for  most  of  the  stratified  tests.  Figure  15 
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shows  a  representative  strip-chart  recording  of  the  internal  wave  disturbance. 
Appendix  B  includes  all  of  the  strip-chart  recordings  obtained. 


Figure  15:  Sample  Internal  Wave  Wire  Strip-Chart  Recording 


In  addition  to  the  computer  and  strip-chart  recordings  of  the  internal 
waves,  VHS  videotape  records  of  the  tests  were  obtained  which  clearly  show  the 
creation  and  propagation  of  the  internal  wave  as  the  ice  floe  moves  through  the 
tank.  Unfortunately,  this  medium  is  not  conducive  for  reproduction  in  a  report. 
However,  the  visualization  of  the  internal  wave  creation  and  behavior  presented 
some  interesting  information.  At  the  speeds  where  the  resistance  “peaked,”  the 
internal  wave  created  by  the  moving  ice  floe  seemed  to  attach  to  the  aft  end  of  the 
floe.  This  is  contrasted  to  the  two  other  conditions:  before  and  after  the  peak.  In 
slower  runs  (before  the  peak),  the  internal  wave  created  was  visible,  but  did  not 
touch  the  ice  floe.  In  faster  runs  (after  the  peak),  the  internal  wave  also  did  not 
touch  the  floe,  and  the  wave  created  diminished  with  increasing  floe  velocity.  See 
Figure  1 6  for  sketches  of  the  floe  with  internal  wave  profiles. 
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Figure  16:  Internal  Wave  Profile  Sketches  for  (a)  Pre-Peak  Resistance, 
(b)  Peak  Resistance,  and  (c)  Post-Peak  Resistance 
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Discussion  and  Preliminary  Analysis 


Non-Stratified  Tests 


Figures  17  and  18  show  the  residual  resistance  coefficients  obtained  from 
the  unstratified  tests  of  Floe  #1  and  Floe  #2  plotted  versus  Reynolds  number.  The 
non-dimensional  Reynolds  number  is  defined  as: 


Re  = 


UL 

V 


where  U  is  the  towing  speed,  L  is  a  characteristic  length  (which  here  is  taken  as 
3  ft  (0.91  m)),  and  v  is  the  fluid’s  kinematic  viscosity. 


Figure  17:  Plot  of  Unstratified  Test  Residual  Resistance,  C,  u^uned» 
vs.  Reynolds  Number  from  Floe  #1 
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Figure  18:  Plot  of  Unstratified  Test  Residual  Resistance,  C  . 
vs.  Reynolds  Number  from  Floe  #2 


One  interesting  observation  of  these  results  is  that  the  residual  resistance 
(which,  in  this  case,  is  directly  proportional  to  form  drag)  is  approximately 
constant  with  Reynolds  number.  As  Bruno  (1992)  discussed,  this  effect  may  be 
explained  by  considering  form  drag  dynamics  of  bluff  bodies.  The  form  drag  of  a 
bluff  body  depends  on  Reynolds  number  only  to  the  extent  to  which  the  boundary 
layer  flow  becomes  turbulent  and  the  flow  separation  pattern  is  complete.  After 
that  point,  there  is  little  dependence  on  Reynolds  number.  Since  the  Reynolds 
numbers  during  the  model  tests  were  sufficiently  high  (order  of  10*),  it  is 
assumed  that  this  critical  point  has  been  surpassed.  Therefore,  the  form  drag,  and 
here  the  residual  resistance,  C„  should  not  vary  significantly  with  Reynolds 
number. 

The  average  residual  resistance  (C,  values  are  plotted  versus  draft 

in  Figure  19.  As  expected,  the  floe  with  the  periodic  bottom  (Floe  #2)  exhibits 
significantly  greater  residual  resistance  than  the  flat-bottom  floe  (Floe  #1)  at  all 
drafts. 
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Stratified  Tests 

The  corresponding  run  number  is  indicated  next  to  each  data  point  for  all 
plots  in  this  section  to  facilitate  easy  reference  to  Table  1  (in  the  Test  Matrix 
section).  The  sketched  curves  in  each  plot  are  intended  to  point  out  noticeable 
trends  in  the  total  resistance  with  respect  to  the  parameters  varied  during  the  tests. 

Figures  20  through  23  present  the  total  resistance  coefficients  from  the 
stratified  tests  (C,  plotted  versus  Reynolds  number  (as  defined  previously). 
An  immediate  observation  of  the  trends  in  the  total  resistance  is  that  the  C,  values 
“peak”  at  certain  Reynolds  numbers  for  each  floe.  Figure  24  shows  the  extracted 
internal  wave  resistance  component,  Cj,,,,  for  Floe  #2.  These  values  were  obtained 
by  subtracting  the  unstratified  test  results  from  the  stratified  test  results.  It  is 
evident  from  this  plot  that  the  internal  wave  effect  diminishes  at  both  low  and 
high  Reynolds  numbers. 

There  are  several  ways  of  presenting  the  data  for  further  interpretation. 
Figures  25  through  28  plot  the  C,  data  versus  a  non-dimensional  floe 
velocity,  U/C,,  where  C,  is  defined  as  the  first-mode  interfacial  phase  speed.  The 
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more  general  C,  relation, 


Q  = 


k  Pi  cothAiAi +p2  coth^2  ’ 
is  the  dispersion  relation  solved  for  a  horizontally  propagating  gravity  wave  in  a 
two-layer,  inviscid  fluid.  In  the  shallow  water,  long  wave  case,  both  kh,  and  kh2 

approach  zero;  cothM  1,2  — ,  and  the  above  equation  may  be  approximated 

^1,2 

by; 

2  gAp  hiht 


C\  = 


P2  hi+hi' 


Non-dimensionalizing  the  results  also  accounted  for  density  variations 
through  the  tests.  Figure  29  shows  selected  data  series  from  each  floe  plotted 
together.  The  similarity  between  the  results  is  intriguing,  considering  the  differing 
bottom  complexities  and  different  T/h,  values  represented. 


Figure  20:  Plot  of  Stratified  Test  Total  Resistance  Coefficient,  C, 


vs.  Reynolds  Number  from  Floe  #1 


37 


lotal  Resistance  Co( 


Figure  21:  Plot  of  Stratified  Test  Total  Resistance  Coefficient,  C, 
vs.  Reynolds  Number  from  Floe  #2 


■  T/h,=  0.46  ♦  T/h,=  0.43  •  T/h,=  0.62  A  T/h,=  0.57 


Figure  22:  Plot  of  Stratified  Test  Total  Resistance  Coefficient, 
vs.  Reynolds  Number  firom  Floe  #3 
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Internal  Wave  Resistance  Coetf,  Ciw 


Reynolds  Number,  \JLN 


Figure  24:  Plot  of  Internal  Wave  Resistance, 
vs.  Reynolds  Number  from  Floe  #2 
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Resistance  Coefficient,  |  Total  Resistance  Coefficient,  C 


Figure  25;  Plot  of  Stratified  Test  Total  Resistance,  C,  sn,tiried» 
vs.  Non-Dimensional  Floe  Velocity,  U/C,  from  Floe  #1 


B 

o 


Figure  26:  Plot  of  Stratified  Test  Total  Resistance,  C, 


vs.  Non-Dimensional  Floe  Velocity,  U/C,  from  Floe  #2 
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sistance  Coeffi 


Figure  27:  Plot  of  Stratified  Test  Total  Resistance,  C, 
vs.  Non-Dimensional  Floe  Velocity,  U/C,  from  Floe  #3 


Figure  28:  Plot  of  Stratified  Test  Total  Resistance,  C,  a„,ified. 


vs.  Non-Dimensional  Floe  Velocity,  U/C,  from  Floe  #4 
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vs.  Non-Dimensional  Floe  Velocity,  U/C,, 
Selected  Data  from  Each  Floe 


Finally,  Figure  30  is  a  plot  of  the  stratified  test  total  resistance  coefficient, 
Q.««iined>  versus  the  draft  to  upper  layer  thickness  ratio,  T/h,.  The  data  presented, 
from  Floe  #1  and  Floe  #3  tests,  illustrates  a  dependence  of  drag  with  respect  to 
floe’s  proximity  to  the  interface  of  the  two  fluids.  Surprisingly,  there  appears  to 
be  a  decline  in  the  total  resistance  as  the  draft  to  upper  layer  ratio  approaches 
unity  for  each  floe.  This  may  be  due  to  modifications  of  the  pressure  disturbance 
below  the  floe  by  the  presence  of  the  fresh  water  /  salt  water  interface. 

The  internal  wave  videotape  recordings  also  supply  some  interesting 
observations.  The  characteristics  of  the  internal  wave  generated  by  the  moving  ice 
floe  appear  less  dependent  on  ice  floe  bottom  topography  than  originally 
expected.  One  of  the  reasons  for  this  occurrence  may  be  that  the  internal  wave  has 
a  much  longer  wavelength  than  the  floe’s  length.  Although  the  pressure 
distribution  immediately  beneath  the  floe  differs  for  each  bottom  topography, 
small  pressure  variations  have  little  effect  on  the  characteristics  of  the  generated 
wave,  and  the  characteristics  of  the  wave  are  dictated  rather  by  the  average 
pressure  disturbance  beneath  the  floe.  It  is  also  interesting  that  the  phenomenon  of 
internal  wave  “attachment”  to  the  aft  end  of  the  ice  floe  (as  illustrated  in 


42 


Figure  30:  Plot  of  Stratified  Test  Total  Resistance  Coefficient,  C, 
vs.  Draft  to  Upper  Layer  Thickness  Ratio,  T/h, 


Figure  1 6,  Test  Results  section)  occurred  with  each  of  the  four  floes  tested,  and 
was  always  coincident  with  the  peak  resistance  coefficient  value.  The  association 
of  maximum  resistance  and  the  “attachment”  phenomenon  is  not  surprising; 
production  of  the  turbulent  mixing  evident  at  the  aft  end  of  the  floe  clearly  adds 
significant  resistance  to  the  moving  floe. 
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Conclusions 


Although  additional  non-stratified  tests  still  need  to  be  completed,  it  may 
be  concluded  now  that  the  contribution  of  internal  wave  creation  on  the  total 
resistance  of  an  ice  floe  moving  in  a  stratified  body  of  water  has  been  proven, 
through  laboratory  experiments,  to  be  significant  and  to  show  definite  trends.  It  is 
apparent  that,  for  a  particular  ice  floe  geometry,  the  additional  resistance  caused 
by  the  wave  creation  “peaks”  at  some  optimal  combination  of  draft  to  upper  layer 
thickness  ratio  and  velocity,  and  then  tapers  off  and  approaches  the  non-stratified 
values. 

The  ramifications  of  the  results  found  in  this  study  on  ice  floe  numerical 
modelling  is  significant.  The  present  methods  for  determination  of  ice-water  drag 
coefficients  should  include  dependencies  on  under-ice  topographies,  as  well  as 
water-body  salinity  profiles.  This  latter  dependence  should  ideally  be  a  dynamic 
one,  as  when  the  ice  floe  travels  from  waters  with  one  salinity  profile  to  another, 
the  corresponding  drag  coefficient  should  also  change. 

It  is  clear  that  the  interrelationship  between  the  factors  relating  to  ice  floe 
movement  is  complex.  Although  certain  correlations  between  ice-water 
resistance  and  several  parameters  have  been  identified,  a  great  deal  of  work  still 
needs  to  be  completed  to  fully  understand  ice  floe  dynamics. 
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Appendix  B 


Internal  Wave  Strip  Chart  Recordings 


Notes  on  internal  wave  wire  strip-chart  recordings: 

•  The  internal  wave  wire  was  mounted  (stationary)  along  the 
side  of  the  tank,  approximately  20  feet  (6  m)  into  the 
data  collection  portion  of  the  test  run. 

•  Each  internal  wave  record  shown  here  is  from  the  period 
during  which  the  ice  floe  model  was  within  the  50  ft 
(15.24  m)  data  collection  length. 

•  The  horizontal  scale  for  all  recordings  is 

1  mm  =  1  sec. 

•  Internal  wave  wire  recordings  were  not  obtained  for 
Runs  23-29  due  to  malfunctions  with  the  instrument. 
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Appendix  C 


Salinity  Profile  Data  and  Plots 


Data  for 

SALINITY  PROFILE  PLOTS 

(Run  #s  8-15) 


(Note:  Salinity  profiles  for  Runs  #1  -7  were  acquired  using  a  point  conductivity  probe. 
Because  of  problems  with  the  probe,  the  data  obtained  from  the  instrument 
was  discarded.) 
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Data  for 

SALINITY  PROFILE  PLOTS 

(Run  #s  16-23) 
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Data  for 

SALINITY  PROFILE  PLOTS 

(Run  #s  24-28,  48-49) 
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SALINITY  PROFILE  PLOTS 

(Run  #s  50-54) 
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